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Popular Summary 

Air pollution has no boundary. Three decades of satellite images provide 
mounting evidence of the intercontinental and hemispheric transport of aerosols, such as 
the unraveling of “Arctic haze”, the sweep of Saharan dust over the Atlantic Ocean, and 
the detection of Asian pollution and dust in North America. Because of poor quality of 
early-generation satellite observations, however, these impressive images only provide 
qualitative or semi-quantitative estimates of long-range transport of aerosols, which is not 
adequate to address important implications of continental scale transport of aerosol for air 
quality, climate change, atmospheric chemistry, and ocean biology applications. 

What can satellite data tell us about the long range transport of pollution beyond 
the pretty pictures? As a result of the much improved measurement accuracy and 
enhanced new capabilities of satellite sensors such as those onboard the A-Train, 
quantitative assessments of aerosol long-range transport from measurements just became 
feasible very recently. Such measurement-based assessments can augment intensive field 
campaigns by expanding temporal and spatial scales and also serve as constraints for 
model simulations. 

In this study we estimated interannual and seasonal variations of mass fluxes of 
pollution aerosol exported from East Asia to the North Pacific basin (East Asia outflow) 
and imported to North America (North America inflow) based on satellite data 
supplemented by in situ measurements. The method takes advantage of MODIS high- 
accuracy measurements of aerosol optical depths over ocean, fine- and coarse-mode 
separately, in combination with satellite measurements of vertical profiles of aerosol and 
water vapor respectively from GLAS and AIRS, and in situ measurements of 
comprehensive aerosol chemistry and physics in the regions. 

Our calculations show that about 18 Tg pollution aerosol is exported from East 
Asia to the northwestern Pacific Ocean each year, of which about 25% reaches the west 
coast of North America. The pollution fluxes are largest in spring and smallest in 
summer, which is mainly controlled by Asian monsoon circulations. During 2002-2005 
the strongest export and import of particulate pollution occurred in 2003, due largely to 
record intense and persistent biomass fires from Siberia in summer and spring. Model 
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simulations of two NASA aerosol models, namely GOCART and GMI, agree quite well 
with the satellite-based assessment of annual pollution fluxes. 

Such measurement-based estimate of trans-Pacific transport is critical to the 
assessment of how pollution from Asia influences climate and air quality in North 
America. It could also play a significant role in formulating effective control strategies in 
North America for combating air quality degradation and climate change. This study is 
the first large-scale measurement-based estimate of intercontinental pollution transport. 
Air pollution is a global problem and the satellite-based approach developed in this study 
can be adapted and extended to study the intercontinental transport of pollution and dust 
in other regions. 
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Abstract: It has been well documented that pollution aerosol and dust from East Asia can 
transport across the North Pacific basin, reaching North America and beyond. Such 
intercontinental transport extends the impact of aerosols for climate change, air quality, 
atmospheric chemistry, and ocean biology from local and regional scales to hemispheric and 
global scales. Long term, measurement-based studies are necessary to adequately assess the 
implications of these wider impacts. A satellite-based assessment can augment intensive field 
campaigns by expanding temporal and spatial scales and also serve as constraints for model 
simulations. Satellite imagers have been providing a wealth of evidence for the intercontinental 
transport of aerosols for more than two decades. Quantitative assessments, however, became 
feasible only recently as a result of the much improved measurement accuracy and enhanced new 
capabilities of satellite sensors. In this study, we generated a 4-year (2002 to 2005) climatology 
of optical depth for pollution aerosol (defined as a mixture of aerosols from urban/industrial 
pollution and biomass burning in this study) over the North Pacific from MODerate resolution 
Imaging Spectro-radiometer (MODIS) observations of fine- and coarse-mode aerosol optical 
depths. The pollution aerosol mass loading and fluxes were then calculated using measurements 
of the dependence of aerosol mass extinction efficiency on relative humidity and of aerosol 
vertical distributions from field campaigns and available satellite observations in the region. We 
estimated that about 18 Tg/year pollution aerosol is exported from East Asia to the northwestern 
Pacific Ocean, of which about 25% reaches the west coast of North America. The pollution 
fluxes are largest in spring and smallest in summer. For the period we have examined the 
strongest export and import of pollution particulates occurred in 2003, due largely to record 
intense Eurasia wildfires in spring and summer. The overall uncertainty of pollution fluxes is 
estimated at about 80%. A reduction of uncertainty can be achieved with a better characterization 
of pollution aerosol through integrating emerging A-Train measurements. Simulations by the 
Goddard Chemistry Aerosol Radiation and Transport (GOCART) and Global Modeling Initiative 
(GMI) models agree quite well with the satellite-based estimates of annual and latitude- 
integrated fluxes, with larger model-satellite differences in latitudinal variations of fluxes. 
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1. Introduction 

1.1. Asian pollution and intercontinental transport 

East Asia is an important source region of a variety of natural and anthropogenic 
aerosols, also called particulate matter (PM). China, the world’s most populous country has been 
undergoing persistent rapid industrialization and urbanization and the burgeoning expansion of 
automobile usage over the last two decades. The nation relies heavily on coal burning to meet 
ever-increasing energy needs {Liu and Diamond, 2005) and consumes large amounts of low 
efficiency biofuels ( Streets and Aunan, 2005). China has more than doubled its emissions in 
recent decades to quickly become one of the largest emitters of aerosols and aerosol precursors 
in the world {Ohara et al., 2007; Richter et al., 2005) and is the largest contributor of sulfur 
dioxide (SO 2 ) to the atmosphere. The Chinese emission of carbonaceous aerosols is nearly the 
sum of North America and Europe emissions {Bond et al., 2004) and the observed concentrations 
for carbonaceous aerosols in mega-cities are 1~3 orders of magnitude higher than those reported 
from developed countries {Wang et al., 2006). While Chinese emissions of SO 2 , carbonaceous 
aerosols and CO showed a decreasing trend during 1996-2000, emissions for these species have 
since been increasing (e.g., Carmichael et al., 2002; Ohara et al., 2007). The increased emission 
of pollutants in East Asia is projected to continue in the coming decades, although an accurate 
projection of the rate of increase is complicated by uncertainties associated with future economic 
growth, air pollution control policy, and bottom-up emission inventory itself, among others 
{Carmichael et al., 2002; Akimoto et al., 2006; Ohara et al., 2007). Dust emissions in the 
northwestern plains of China are comparable with that of the Saharan dust. Dust storms 
frequently blanket a large area of northern China in the spring season and have become 
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increasingly intense over the past few years after a decreasing trend from 1969-1999 ( Zhang et 
al., 2003). 

Such large emissions of aerosols may have brought about substantial impacts to China’s 
climate (Huang et al., 2006), hydrological cycle ( Menon et al., 2002), solar radiation reaching 
the surface (Qian et al., 2007), crop yields (Chameides et al., 1999), and human health (Xu et al., 
1998). Furthermore, Asian aerosol could impose far-reaching environmental impacts at regional, 
hemispherical and even global scales, because of mounting evidence for intercontinental and 
even hemispheric transport provided by long-term surface monitoring networks, in-situ 
measurements from intensive field campaigns, and satellite observations backed by model 
simulations. Dust storms and industrial emissions have long been observed to influence the 
North Pacific (Duce et al., 1980), to contribute to particulate matter level on the western 
boundary of North America (Andreae et al., 1988; Jajfe et al., 1999; Husar et al., 2001 ; 
Wilkening et al., 2000; McKendry et al., 2001; Jajfe et al., 2005) and on the East Coast after 
traversing the continent (Biscaye et al., 2000; VanCuren and Cahill, 2002), and to even affect 
Europe (Stohl et al., 2007; Grousset et al., 2003). Model simulations also suggest substantial 
trans-Pacific transport of aerosols (Chin et al., 2004, 2007; Heald et al., 2006; Hadley et al., 
2007). 


East Asia is one of the main warm conveyor belts (WCB) inflow regions, where 
airstreams transport the boundary layer air into the upper troposphere on time scales of 1-2 days 
(Stohl, 2001; Eckhardt et al., 2004). Other vertical lifting mechanisms include frontal and 
postfrontal convection, orographic lifting, and the atmospheric boundary layer (ABL) turbulent 
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mixing. As a result, emissions from Asia can experience rapid vertical transport ( Stohl et al., 
2002). Pollutants lifted to the free troposphere travel initially poleward and eastward by mid- 
latitude storm tracks, turn towards the equator thereafter, and then reach mainly the middle and 
upper troposphere of the northeastern Pacific. Such high-level transpacific transport takes about 
a week {Holier et al, 2005). The polluted Asian air over the Northeastern Pacific can then be 
transported downward to influence the North American ABL through subsidence associated with 
subtropical Pacific highs fHeald et al., 2003) and through meso-scale mountain circulations in 
Western North America {McKendry et al., 2001). Air pollutants can also be transported across 
the North Pacific at low altitudes through ABL outflow associated with cold fronts {Liang et al., 
2004; Liu et al., 2003), with transit times 1-2 weeks longer than the high-level trans-Pacific 
transport {Holzer et al., 2005). Specific mechanisms responsible for trans-Pacific transport 
depend strongly on season and location {Eckhardt et al., 2004; Liang et al., 2004; Holzer and 
Hall, 2007; Dickerson et al., 2007). This in combination with variations in pollutant emissions 
could result in pronounced seasonal and interannual variations of pollution outflow and inflow. 

The regional and global impacts of Asian aerosols are largely unknown but the 
implications are quite powerful. A tracer simulation suggests that Asian CO emissions with an 
assumed atmospheric lifetime of 5 days (comparable to lifetimes of tropospheric aerosol) could 
contribute to the total CO column over North America by about 30% {Stohl et al., 2002). 

Because particulates with a diameter less than 2.5 pm (PM2.5), can efficiently penetrate into the 
lungs and cause cardiovascular and respiratory diseases, increasingly stringent standards for 
PM2.5 have been enforced by the U.S. Environmental Protection Agency (EPA). While most 
pollution sources are local, the trans-boundary and intercontinental transport of PM2.5 from Asia 
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may compromise efforts of some western U.S. states to attain air quality goals through domestic 
and regional emission controls (e.g., Heald et al., 2006; Chin et al., 2007; Fairlie et al., 2007). 
Aerosols can perturb the radiative budget of the earth-atmosphere system directly through 
scattering and absorbing radiation and indirectly through altering cloud formation and cloud 
properties. The deposition of soot in snow and ice can also reduce surface reflectance and hence 
perturb the surface energy budget ( Hansen and Nazarenko, 2004). Linked through these 
processes, aerosols can thereby affect atmospheric circulations and precipitation. It has been 
suggested that an increase of particulate air pollution may reduce orographic precipitation in the 
mountain ranges of the western United States ( Rosenfeld and Givati, 2006; Jirak and Cotton, 

2006) , causing a significant shortage of water resources in the western states. Asian dust particles 
have been observed to affect the formation of clouds over the western US ( Sassen , 2002). A 
satellite-detected intensification of the North Pacific mid-latitude storm tracks has been 
associated with an increase of aerosol emission in Asia in the last two decades {Zhang et al., 

2007) . This suggests a possible feedback influence of long-range transport of aerosols, because 
the North Pacific storm tracks effectively transport pollution and dust from East Asia to North 
America. 

1.2. A satellite perspective 

A better understanding of aerosol transport, transformation, and distribution is required to 
assess a variety of aerosol impacts on a wide range of scales and formulate an effective strategy 
for combating air pollution in the U.S. Due to the large emissions and rate of increase in Asia, 
the potential impacts on the environment are relatively easy to detect and can be studied over a 
large dynamic range. As such, in the past two decades, the aerosol related scientific and 
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environmental issues in the North Pacific basin have been a focus of several large scale 
international experiments, including the Western Pacific Exploratory Missions (PEM-West A 
and B) ( Hoell et al., 1996, 1997), Aerosol Characterization Experiment- Asia (ACE- Asia) 

( Huebert et al., 2003), TRAnsport and Chemical Evolution over the Pacific (TRACE-P) {Jacob 
et al., 2003), Photochemical Ozone Budget of the Eastern North Pacific Atmosphere (PHOBEA) 
( Bertschi , et al., 2004), Intercontinental Transport and Chemical Transformation (ITCT 2K2), 
and Pacific Exploration of Asian Continental Emission (PEACE) {Parrish et al., 2004), 
Intercontinental Chemical Transport Experiment - Phase B (INTEX-B), among others. During 
some of these intensive field campaigns, aerosols have been studied in great detail using in-situ 
and remote sensing observations of physical and chemical properties from multiple platforms 
such as aircraft, ship, and ground networks and model simulations as well. These missions have 
been providing comprehensive snapshots of regional aerosols that can be intercompared to help 
understand the complex interactions of aerosols within atmospheric systems. 

For all of their advantages, field campaigns are inherently limited by their relatively short 
duration and small spatial coverage. Therefore they are not adequate to assess the temporal and 
spatial variations of trans-Pacific transport of aerosols. Satellite remote sensing of aerosols can 
augment these missions by expanding temporal and spatial scales to generate measurement based 
estimates of aerosol intercontinental transport. In addition, these estimates can serve as 
constraints for model simulations. Satellite imagers not designed specifically for aerosol research 
have been providing a wealth of evidence for intercontinental and hemispheric transport of 
aerosols in a qualitative or semi-quantitative manner for nearly three decades (e.g., Lyons et al., 
1978; Chung, 1986; Herman et al., 1997). Quantitative estimates of aerosol mass and transport 
were also attempted {Fraser, 1976; Fraser et al., 1984). However these early estimates were 
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subject to large uncertainties because of poor accuracy of early satellite measurements. 
Quantitative assessments of aerosol intercontinental transport became feasible only recently as a 
result of the much improved measurement accuracy and enhanced new capabilities of satellite 
sensors (King et al., 1999; Kaufman et al., 2002), as well as a constellation of multiple satellites 
(such as “A-Train”) with complementary measurement capabilities. Although current and near- 
future satellite sensors can not directly measure aerosol composition, measurements can be used 
to categorize aerosol types in terms of particle size (i.e., fine vs coarse), shape (i.e., spherical vs 
non-spherical), and absorption (i.e., scattering vs absorbing) because of implementation of multi- 
wavelength, multi-angle, and polarization measurements (Higurashi and Nakajima, 2002; Tanre 
et al., 2001). Algorithms have been developed to quantitatively estimate the optical depth of 
pollution aerosol and dust (e g., Kaufman et al., 2005a, b). The new capabilities associated with 
these passive sensors are being further enhanced by emerging lidar measurements in recent years 
that are providing essential information of aerosol vertical distributions ( Spinhime et al., 2005; 
Winker et al., 2003). It is currently feasible to develop a satellite-based approach supplemented 
by in situ measurements and model simulations to quantify aerosol intercontinental transport. 

In this study we make our first attempt at this approach by estimating the mass fluxes of 
pollution aerosol exported from East Asia to the North Pacific basin and imported to North 
America based on satellite data. For this purpose, the satellite observed aerosol is categorized 
into three generic types based on its sources, i.e., maritime aerosol, mineral dust, and pollution 
aerosol. Here pollution aerosol is considered to be a mixture of aerosols from urban/industrial 
pollution and biomass burning, which is predominately in sub-micron size or fine-mode. All 
biomass burning is assumed to be man made here. The method takes advantage of the Moderate 
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resolution Imaging Spectroradiometer (MODIS) high-accuracy measurements of aerosol optical 
depths over ocean, fine- and coarse-mode separately, in combination with satellite measurements 
of aerosol vertical profiles from the Geoscience Laser Altimeter System (GLAS), and in situ 
measurements of comprehensive aerosol chemistry and physics from intensive field campaigns 
in the regions. The estimated aerosol import and export fluxes are compared with the Goddard 
Chemistry Aerosol Radiation Transport (GOCART) and Global Modeling Initiative (GMI) 
simulations. The satellite observations of aerosol and aerosol models are described in Section 2. 
In Section 3, an approach to deriving pollution fluxes from satellites is presented. Seasonal and 
interannual variations of the satellite-estimated export and import of pollution aerosol are 
discussed and compared with the model simulations in Section 4. Uncertainties associated with 
the flux estimate and potential of reducing uncertainties by integrating the Cloud- Aerosol Lidar 
with Orthogonal Polarization (CALIOP) and other A-Train observations are also discussed. 
Major results are summarized in Section 5. 

2. Description of Aerosol Data and Models 

2.1. MODIS aerosol optical depths 

MODIS aboard both Terra and Aqua is making near global daily observations of 
atmospheric aerosols since February 2000 (Terra) and July 2002 (Aqua). MODIS uses seven 
wavelength channels (between 0.47 and 2.13 pm) to retrieve aerosol properties over cloud and 
surface-screened areas. Separate algorithms are implemented over land and ocean ( Kaufman et 
al., 1997; Tanre et al., 1997; Renter et al., 2005; Levy et al., 2007). MODIS products are widely 
used by the community for a variety of purposes, including descriptions of the regional, seasonal 
and global distribution of aerosols and their relationship to other pollutants (Yu et al., 2003; Chin 
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et al., 2004), studies of the effect of aerosol on atmospheric chemistry and local air pollution 
(Engel-Cox et al., 2004; Al-Saadi et al., 2005), measurements of the aerosol radiative forcing of 
climate (Yu et al., 2004; Renter and Kaufman, 2006), studies of the aerosol interaction with the 
meteorological field and with clouds (Koren et al., 2004; Kaufman et al., 2005d), and synergy 
with other sensors to enhance standard MODIS retrievals and offer new standalone products 
(Kaufman et al., 2003a,b; Leon et al., 2003). 

Because of its wide spectral range over ocean, MODIS has the unique capability to 
retrieve not only spectral aerosol optical depth (x) at seven wavelengths from 0.47 to 2.13 (am 
with great accuracy, i.e., ±0.03±0.05x (Remer et al., 2002, 2005), but also the quantitative 
aerosol size parameters (e.g., effective radius and fine-mode fraction) (Kaufman et al., 2002; 
Remer et al., 2005; Kleidman et al., 2005). The fine-mode fraction (FMF) is a measure of the 
contribution of fine-mode aerosols to the optical depth. Because anthropogenic aerosols are 
predominately fine-mode or in submicron range, FMF can be used as a tool for separating 
anthropogenic aerosol from dust, as demonstrated in some previous studies. Kaufman et al. 
(2005a) developed an algorithm that uses MODIS/Terra aerosol optical depth (AOD) and FMF 
to distinguish dust from pollution and maritime aerosols and evaluate the column concentration, 
transport and deposition of dust over the Atlantic Ocean. The results confirm seasonal variations 
of magnitude and location of African dust transport across the Atlantic Ocean. In particular, the 
estimated deposition of 50 Tg dust in the Amazon, much larger than previous estimates, may 
explain the paradox between the need of nutrition by the Amazon forest and the source of the 
nutrition. A similar algorithm has been employed to estimate anthropogenic and natural 
components of aerosols (Kaufman et al., 2005b), which has since inspired the community to 
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further explore the use of satellites to quantify aerosol radiative forcing by anthropogenic aerosol 
(e.g., Anderson et al., 2005a, Bellouin et al., 2005; Christopher et al., 2006; Yu et al., 2006). 

2.2. ICESat and CALIPSO aerosol profiles 

GLAS was launched in early 2003 aboard ICESat to become the first polar orbiting 
satellite lidar with the ranging capability to serve as both a precision surface elevation altimeter 
and an atmospheric lidar ( Spinhirne et al., 2005). CALIPSO is a joint U.S. and French satellite 
mission with an expected 3 year lifetime ( http://www-calipso.larc.nasa.gov/ ). Launched on April 
28, 2006, CALIPSO combines an active lidar instrument (CALIOP) with passive infrared and 
visible imagers ( Winker et al., 2003). GLAS and CALIOP both use 532 nm and 1064 nm 
wavelengths to probe the vertical structure and properties of thin clouds and aerosols with a 
resolution of ~60 m vertically and ~10 km along track. The 532 nm wavelength of CALIOP also 
has a capability of detecting polarization, which can facilitate the separation of non-spherical 
dust from spherical particles. These satellite sensors cover the globe between ~80°N and ~80°S 
latitudes, with an orbital repeating cycle of about 16 days. Although they are nadir track only 
sensors, useful global distributions can emerge for monthly compilations involving hundreds of 
orbits. 


Since February 2003, GLAS has operated during specific time segments of about 4 to 8 
weeks duration due to the limited laser lifetime. Although limited in time coverage, the GLAS 
data provide, for the first time, the vertical distributions of aerosols backscatter at 532 nm up to a 
40 km height on a global scale ( Spinhirne et al., 2005). While similar products will be derived 
from CALIPSO measurements that are expected to be continuous in temporal coverage, the full 
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calibrations and aerosol-cloud separation have not yet been implemented at the time of 
conducting this study. In this study, we will examine GLAS data only and discuss the potential 
of CALIPSO data for better determining the major transport heights of pollution plumes and a 
reduction of uncertainty associated with flux estimates in Section 4.3. 

2.3. GOCART and GMI aerosol simulations 

The global model GOCART is driven by assimilated meteorological fields from the 
Goddard Earth Observing System Data Assimilation System (GEOS DAS) ( Chin et al., 2002). 
GOCART simulates the major aerosol types: sulfate, dust, black carbon, organic carbon, and sea- 
salt. Emissions from anthropogenic, biomass burning, biogenic, and volcanic sources and wind- 
blown dust and sea-salt are included in the model. Processes represented in GOCART are 
chemistry, convection, advection, boundary layer mixing, dry and wet deposition, gravitational 
settling, and hygroscopic growth of aerosol particles. Details of GOCART and evaluation of its 
results against observations are documented in a number of publications (e.g., Chin et al., 2000a, 
2000b, 2002, 2004; Ginoux et al., 2001, 2004). 

The Global Modeling Initiative ( GMI) was established to understand and reduce 
uncertainties in assessing the impact of various natural and anthropogenic perturbations on 
atmospheric composition and chemistry. The GMI aerosol module, adapted from the University 
of Michigan/LLNL IMPACT model ( Liu et al., 2007), is similar to the GOCART model in such 
aspects as fossil fuel emissions of SO 2 , BC, and OC, and the advection core. On the other hand, 
the two models differ in volcanic emissions, aqueous phase reactions, vertical diffusion, and dry 
deposition and wet scavenging schemes, among others. As such GOCART and GMI simulated 
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mass and flux of pollution aerosol can be different both in the horizontal and vertical, though 
both models are driven by the same meteorological fields from GEOS-4 in this study. 

For a comparison with MODIS-based estimates of pollution flux, we use a sum of model 
simulated black carbon, organic matters, and sulfate as a proxy for pollution aerosols. While this 
sum is higher than the modeled “pollution” aerosols since there are natural components of sulfate 
and organic matter, the bias could be partially offset by an exclusion of nitrate and ammonium in 
the models that is largely anthropogenic. Note also that for the MODIS estimate we have 
assumed all biomass burning smoke is man-made pollution, which is somewhat consistent with 
the pollution proxy assumption made for the models. 

3. Approaches to estimating pollution fluxes from satellites 

In this study we extend the method and analysis in Kaufman et al. (2005a, b) to estimate 
fluxes of pollution aerosol exported from East Asia and imported to North America. A three-step 
method integrates datasets from multiple satellite sensors and in situ measurements from field 
campaigns to generate a four-year (2002-2005) climatology of trans-Pacific pollution aerosol 
transport, as discussed in the following. 

3.1. Deriving pollution aerosol optical depth from MODIS retrievals 

We use the method developed by Kaufman et al. (2005a,b) to derive pollution aerosol 
optical depth (x p ) according to: 

_ _ (f ~ f ~ (f m ~ f d )^m / 1 \ 

(/,-/,) 0) 
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where % and f represents aerosol optical depth and fine-mode fraction respectively, both at 550 
nm. Subscripts p, d, and m denote pollution, dust, and maritime aerosol, respectively. Here we 
assume that volcanic aerosols are negligible. Such assumptions may bias our estimate of 
pollution optical depth to a high value in years following major volcano eruptions. We applied 
the method to Collection 4, Level 3 daily MODIS/Terra aerosol data (at l°xl° resolution) so that 
empirical coefficients in (1), f m , f d , f p , and T m derived in Kaufman et al. (2005a, b) are valid. 

We modify the method developed in Kaufman et al. (2005a, b) by taking account of 
seasonal and spatial variations of fine-mode fraction for maritime aerosol (f m ). Over remote 
oceans, boundary layer aerosols consist of sea salt particle generated from bursting bubbles, 
sulfate from oxidation of dimethyl sulfide (DMS), and organic matter. While DMS-oxidized 
sulfates and organic particles are in the sub-micron size range, the sea salt aerosols have much 
broader size distributions {Bates et al., 2001) and the fine-mode sea salt is a significant 
contributor to sea salt optical depth and also an important component of fine-mode maritime 
aerosol optical depth. These marine-generated fine-mode aerosols should depend on a variety of 
atmospheric and oceanic parameters, such as ocean color, sea-surface temperature (SST), near- 
surface wind speed, and atmospheric oxidizing capacity. This complexity gives rise to large 
seasonal and geographical variations in the maritime FMF (ranging from 0.2 to 0.8), as 
evidenced consistently in a MODIS-derived proxy for maritime aerosol (e.g., AOD less than 0.1) 
and GOCART and GMI simulations of purely maritime aerosol (Yu et al.. Seasonal and 
geographical variations of maritime aerosol fine-mode fraction and implication for deriving 
MODIS anthropogenic aerosol, in preparation). A use of the seasonally and geographically 
varying maritime FMF decreases the pollution aerosol optical depth by as much as 30% in 
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coastal areas of the North Pacific and in summer, in comparison to that using the constant f m 
(Kaufman et al., 2005 a,b). 

3.2. Determining major transport heights of pollution aerosol 

To transform aerosol optical depth to aerosol mass loading and flux, major transport 
heights of pollution aerosol are needed to determine the relative humidity and wind speed 
representative of satellite observed pollution plumes. Passive sensors like MODIS acquire 
columnar aerosol properties which give essentially no information on vertical distributions of 
aerosol. In this study we determine the transport heights of pollution aerosol based on aircraft 
measurements in the region, available surface and satellite lidar observations, and model 
simulations. 

Observations over the northwestern Pacific during ACE- Asia, TRACE-P, and PEM-West 
B indicate multiple pollution layers from the surface up to the middle troposphere with a 
significant fraction of pollution aerosol mass above the ABL (e.g., Bahreini et al., 2003). For 
example, S 04 = and NH/ profiles over Northwestern Pacific indicate that about one third of 
sulfate mass is above 2 km (Dibb et al., 2003; Chin et al., 2003). Carbonaceous aerosols were 
even more concentrated than sulfate in the upper troposphere ( Huebert et al., 2004). A global 3- 
D chemical transport model showed that the trans-Pacific transport of Asian anthropogenic 
aerosols takes place mainly in 900-700 hpa layer because of scavenging during transport either in 
the boundary layer or during lifting to the upper troposphere (He aid et al., 2006). Measurements 
over the Northeastern Pacific Ocean and Western North America indicate persistent Asian 
continental influences in the lower free troposphere (VanCuren et al., 2005). While Siberian 
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boreal fire emissions and Asian dusts tend to elevate the aerosol extinction mainly above the 
ABL (from 2-3 km up to 6 km) over the Northeastern Pacific, the enhancement of aerosol 
extinction in the ABL results largely from transport of industrial pollution aerosol ( Bertschi et 
al., 2004; Price et al, 2003; de Gouw et al., 2004). Air masses from Southeast Asia and China 
are generally observed at higher altitudes than air from Japan and Korea ( de Gouw et al., 2004). 
Tracer model simulations suggest that Asian pollution transported to the west coast of the 
Northwestern U.S. was largely confined to the lowest 3 km layer, whereas the isentropic 
trajectory indicates a somewhat higher transport height ( Jaffe et al., 1999). 

Multi-year measurements from the ground-based micro-pulse lidars in the Korean 
peninsula {Kim et al., 2007) indicate that, on an annual basis, roughly 2/3 of columnar AOD is 
contributed by aerosols below 2 km. In comparison to summer and fall measurements, aerosol 
extinction in spring and winter is enhanced above the boundary layer and below about 6 km, 
which is presumably attributable to seasonal variations of Asian dust storm s . Lidar 
measurements over Japan and Korea during the ACE- Asia give the similar vertical distribution 
for polluted cases (Quinn et al., 2004). Observations of aerosol profiles from GLAS are also 
examined for two available periods. Figure 1 shows profiles of the frequency and average 
backscatter of GLAS detected aerosol layers during September-November, 2003 (Figure lb) and 
February-March, 2004 (Figure lc) in four coastal regions as illustrated in Figure la. While 
aerosol layers are largely detected in the boundary layer in both periods and all regions, some 
enhanced layers with backscatter comparable to that in the ABL are indeed detected in the free 
atmosphere. Such enhancement in the free atmosphere is more pronounced in February-March 
than in September-November, presumably because of more frequent dust plumes in late winter 
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and spring that generally transports at higher altitudes than pollution plumes. Note also that the 
detected aerosol in the ABL could also be partially attributable to the maritime aerosol. It is thus 
difficult to characterize pollution aerosol from only GLAS aerosol measurements. 


To be consistent with the observations and models discussed above, we assume that the 
aerosol optical depth in each of three layers, namely 0-1 km, 1-2 km, and above 2 km (denoted 
respectively as LI, L2, and L3), accounts for one third of the columnar pollution AOD. That is, 
67% of pollution AOD is assume to reside in the lowest 2 km layer (approximately, ABL) and 
33% above the ABL. This is a simplified representation for the complicated layering structure of 
aerosols in reality. Seasonal, interannual, and geographical variations of pollution transport 
heights are not considered in this study because of the paucity of available observations. 
Therefore seasonal, interannual, and geographical variations of pollution fluxes discussed in this 
study result from a combination of the MODIS observed variations of column aerosol mass 
loading and variations of wind and humidity profiles. 


3.3. Calculating pollution aerosol mass loading and flux 

Pollution aerosol optical depth (x p ) derived above is for hydrated pollution aerosol in 
ambient conditions measured by MODIS. In each layer (1), x p can be converted to dry mass 
loading M p as follows: 


M P {1) 


T P V) 

f[RH(l)]MEE p 


( 2 ) 


where MEE P (m 2 g' 1 ) is dry mass extinction efficiency of pollution aerosol; and f(RH) is a 
function that accounts for an increase of aerosol extinction with increasing relative humidity 
(RH), depending on chemical composition with monotonic (smoothly varying) or deliquescent 
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(step change) growth ( Quinn et al., 2005; Carrico et al., 2003). For MEE P we use a value of 4 
m g' based on measurements during ACE- Asia campaign for continental outflows with weak 
dust influence and at low RH ( Bates et al., 2006 and references therein). For f(RH) we use an 
empirical function that is an average of the monotonic and deliquescent growth curves observed 
during ACE- Asia for pollution aerosol ( Carrico et al., 2003). Measurements of RH from the 
Atmospheric Infrared Sounder (AIRS) ( Aumann et al., 2003) are used as input to quantify the 
hygroscopic growth of aerosols and convert the measured optical depth in ambient conditions to 
the (dry) mass loading. To comply with vertical resolution of AIRS datasets, we use RH at 925, 
850, and 700 hPa to derive pollution mass loading in LI, L2, and L3 layer, respectively. 

The daily mass loading (gm" 2 ) for pollution aerosol derived at individual l°xl° grids is 
aggregated to the monthly average (M p ) over the flux calculation boxes. Then monthly average 
zonal winds at 925, 850, and 700 hPa from the Goddard Earth Observing System Data 
Assimilation System (GEOS DAS) with current version of GEOS-4 are used to calculate 
pollution fluxes (gs' 1 ) in LI, L2, and L3 layer, respectively. For a segment with a length of L (m) 
that is parallel to longitudes, we use east-west wind component U to calculate the columnar 
fluxes as follows: 

F„=±M r (im)*L ( 3 ) 

/=1 

The calculated monthly fluxes are aggregated into seasons and 10° latitude sections for analysis 
in the next section. 

4. Results and Discussion 


17 



Transpacific Transport of Pollution Aerosol from Satellite 


Yu et al. 


We applied the method discussed in Section 3 to MODIS/Terra Level 3 daily l°xl° data from 
2002-2005 to estimate fluxes of pollution aerosol exported from East Asia to West Pacific (WP) 
and imported to the West Coast of North America over East Pacific (EP). The flux estimation 
was done across meridional planes with a width of 10° in longitude, as illustrated in Figure 2. 
Over West Pacific, the meridional plane is centered at 130°E south of 40°N and at 140°E north 
of 40°N. On the other side of Pacific, the meridional plane is centered at 130°W. We integrate 
the calculated flux from 20°N to 60°N over WP to represent the East Asia outflow. Given that 
the transpacific transport usually shifts poleward after leaving East Asia and the transport in 
subtropical (20°-30°N) EP is predominately controlled by easterlies, fluxes are integrated from 
30°N to 60°N over EP to represent North America inflow. In the following, seasonal and 
latitudinal variations of MODIS -estimated export and import fluxes are discussed and compared 
with GOCART and GMI simulations in Section 4.1. Interannual variability of pollution fluxes 
are examined for 2002-2005 period in Section 4.2. In Section 4.3, uncertainties associated with 
flux calculations are estimated and a potential of reducing the uncertainty is discussed in terms of 
integrating A-Train observations. 

4.1. Seasonal and latitudinal variations of pollution fluxes 

Since most field experiments and promoted modeling studies of aerosols have been 
focusing on late winter and spring (e.g., Tan et al., 2002; Chin et al., 2004; Heald et al., 2006; 
Hadley et al., 2007), seasonal variations of transpacific transport of pollution aerosol are not well 
understood. Figure 3 shows the MODIS -estimated seasonal pollution fluxes for 2004. Clearly, 
both East Asia outflow and North America inflow peaked in spring, with the flux of 6.8 and 1 .7 
Tg/season, respectively. The outflow and inflow are weakest but not negligible in summer, with 
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a magnitude of about 1/3 of the springtime maximum and 10% of the annual flux. Winter and 
autumn appear to be transitional seasons when pollution fluxes are in between the springtime 
maximum and the summertime minimum. Over the western Pacific, the wintertime pollution flux 
is about 30% higher than that in autumn. Such seasonal variation for pollution aerosol fluxes is 
similar to simulations of carbon monoxide ( Liu et al., 2003). 

The observed seasonal variations of pollution aerosol fluxes are determined by a 
combination of meteorological conditions, emissions, chemistry and removal processes. In East 
Asia, emissions of pollution aerosols and their precursors have some seasonal variations. In 
winter the use of coal in northern China is more than those in other seasons because of heating 
demands in winter. On the other hand, photochemical production of sulfate aerosol from the 
combustion generated SO 2 is slower in the winter because of lower oxidant levels and slower 
reaction rates. Biomass burning shows strong seasonal variations; In South and Southeast Asia, 
the burning occurs mainly in spring, while in high latitude Eurasian region, the burning peaks in 
spring and summer (e.g., Giglio et al., 2006; Bian et al., 2007). Presumably the Asian monsoon 
system would play a much larger role in regulating the seasonality of trans-Pacific transport, as 
suggested by a number of previous studies. The WCBs are stronger and occur more frequently in 
winter and spring ( Eckhardt et al., 2004) and the springtime dry convection can also be an 
important lifting mechanism over China ( Dickerson et al., 2007). This in combination with 
strongest mid-latitude westerlies results in the strongest pollutant transport in spring. The 
weakest pollutant transport in summer is a result of an offset of the westward transport by the 
eastward transport, particularly south of 30°N and the largest aerosol removal associated with 
summer monsoon circulation ( Holzer et al., 2005). 
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On an annual basis, the MODIS -estimated pollution outflow and inflow is 16.0 and 4.2 
Tg/year, respectively (Table 1). This suggests that about a quarter of East Asia outflow of 
pollution aerosol can reach the west coast of North America, if there is no additional contribution 
of pollution aerosol from elsewhere during trans-Pacific transport. While this would be a 
reasonable assumption on an annual basis, the assumption may overestimate Asia’s contribution 
to North America in summer season and certain years because of potential contributions of 
biomass burning smoke from Alaska. The smoke from Mexico and Central America in spring is 
unlikely to contribute significantly to the calculated North America inflow over 30°-60°N. By 
comparison, GOCART and GMI simulations give an outflow of about 10% smaller than the 
MODIS estimate. On the other hand, both GOCART and GMI give 4.8 Tg/year for the North 
America inflow, i.e., 15% larger than the MODIS-based estimate. The model simulations suggest 
that about one third of East Asia outflow can reach the west coast of North America, which is 
somewhat larger than the MODIS-based estimate. Figure 4 shows that on a seasonal basis 
satellite estimated and model simulated outflow and inflow fluxes are still in reasonably good 
agreement. The largest satellite-model differences occurred in spring for the East Asia outflow 
(as large as about 30% or 2 Tg) and in summer for the North America inflow (as large as 70% or 
0.5 Tg). 

Meridional variations of annual pollution aerosol fluxes are compared between satellite 
estimate and model simulations, as shown in Figure 5. For East Asia outflow, MODIS estimated 
the maximum flux of 7.6 Tg/year in the 30°-40°N segment, followed by 5.7 Tg/year in the 40° - 
50°N segment. The flux integrated over the 30°-50°N segment accounts for about 80% of the 
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total East Asia outflow. A remaining 20% of the outflow is nearly evenly distributed in the 20°- 
30°N and 50°-60°N segments. For comparison, GOCART simulation seems to shift the pollution 
outflow northward. The maximum outflow occurs in the 40°-50°N segment, larger than the 
corresponding MODIS estimate by about 20%. On the other hand, the GOCART simulation of 
outflow is larger than the MODIS estimate by 75% in the 50°-60°N segment but smaller by 60% 
in the 30°-40°N segment. GMI simulations give that major outflow path is situated in the 30°- 
50°N segment, consistent with the MODIS estimate. However the magnitude of flux is about 
25% lower than the MODIS estimates. On the other hand, the GMI flux of 2.9 Tg/year in the 
20°-30°N segment is a factor of more than 2 larger than the MODIS estimate. 

For the North America inflow, MODIS estimates the highest flux of 1.8 Tg/year in the 
40°-50°N segment, followed by 1.4 Tg/year in the 50°-60°N segment and 1.1 Tg/year in the 30° - 
40°N segment. This suggests a general poleward shift of pollution plume during the trans-Pacific 
transport. In the subtropical segment (20°-30°N), the pollution flux is negative, suggesting a 
westward transport of pollution from North America. GOCART fluxes are higher in the 40°- 
60°N segments but lower in the 20°-40°N segments than the MODIS estimates. Similar to that 
over West Pacific, GMI seems to shift the inflow toward lower latitudes, with its flux in the 30°- 
40°N segment substantially larger than both MODIS and GOCART fluxes. In the 20°-30°N 
subtropical segment the GMI flux has a small but positive value, suggesting a weak eastward 
transport. This difference is likely to result from different transport heights in conjunction with 
the shift of westerly boundary layer winds to easterlies in the free troposphere. 
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The above analysis suggests that in 2004 about 16Tg pollution aerosol was exported from 
East Asia to the northwestern Pacific Ocean through mainly the 30°-50°N segment. The 
pollution flux was largest in spring and smallest in summer. The summer flux was about 30% of 
the springtime flux. About one quarter of the pollution export from Asia may have reached the 
west coast of North America in that year. This substantial trans-Pacific transport of pollution is 
arguably mostly originating from Asia, with some contributions from Europe, Africa, and 
elsewhere depending on season and location {Newell and Evans, 2000; Bey et al., 2001). The 
satellite-based estimates agree with GOCART and GMI simulations remarkably well over broad 
meridional segments and on an annual basis. The comparisons also show reasonably good 
agreements over 10° meridional segments and on a seasonal basis. In comparison to MODIS 
estimates, GOCART simulations generally shift the pollution transport northward while GMI 
simulations do otherwise. 

4.2. Interannual variability of pollution fluxes 

While some interannual variations or trends of industrial emissions in East Asia have 
been suggested (e.g., Carmichael et al., 2002; Richter et al., 2006), biomass burning emissions 
that are controlled by meteorological and biospheric conditions as well as human activities are 
likely to show much larger interannual variations {Duncan et al., 2003). On the other hand, the 
interannual variability of Asian outflow resulting from large-scale variations in atmospheric 
circulation is the lowest, in comparison to other regions {Liu et al., 2005). It is thus expected that 
pollution fluxes exported from East Asia and imported to North America would show some 
interannual variations resulting mainly from variations in biomass burning smoke from Eurasia. 
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We have applied the method to MODIS/Terra data for other years and derived the four- 
year climatology (2002-2005) of pollution fluxes exported to the West Pacific and imported to 
the East Pacific. Figure 6 shows interannual variations of pollution aerosol optical depth in 
MAM and JJA. In MAM and JJA, 2003, pollution aerosol optical depth in the high-latitude 
North Pacific (40°-60°N) was elevated to a level that is more than a factor of 2 larger than in 
other years. This is corroborated with elevated TOMS absorbing aerosol index that suggests 
intense biomass burning smoke (not shown). In 2003, biomass fires burned a large area of south 
and east Russia, starting from May through the summer ( Goldammer et al., 2004; Wotawa et al., 
2006). Both the fire pixel count and burned area are a factor of 3-4 larger than that in 2004 
( Wotawa et al., 2006). In fact the burned area in 2003 was the largest in a decade. The pollution 
aerosol optical depth in high-latitudes (40°-60°N) is about a factor of 3 larger in 2003 than 2004. 
On the other hand, boreal forest fires in Alaska and Canada in summer were more intense in 
2004 than in 2003. The burned area and fire counts were a factor of 2-3 larger in 2004 than in 
2003. Although such boreal fires mainly influence the North America continent in the eastern 
and southeastern U.S. ( Wotawa and Trainer, 2000), there is a discemable influence on the 
surrounding northeastern Pacific, as shown in the MODIS pollution aerosol map. In the mid- 
latitude West Pacific, pollution largely from China influenced the region to the greatest extent in 
MAM and JJA, 2005. In the tropical East Pacific, the spring-time biomass burning smoke from 
Central America was weakest in 2004. 

Table 2 compares annual pollution fluxes of East Asia outflow and North America 
inflow from 2002 to 2005. The annual pollution flux varies by as much as 25% for the outflow 
and 45% for the inflow. Both the outflow and inflow were strongest in 2003, because of record 
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intense biomass burning emissions from spring to summer in Eurasia. About 20 Tg pollution 
aerosol was transported from East Asia to the northwestern Pacific, of which 28% reached the 
west coast of North America. This is followed by 2005 when 18.6 Tg pollution aerosol was 
exported from East Asia and 4.1 Tg was imported to North America. As discussed earlier, the 
pollution AOD in 2005 was higher than that in other years south of 40°N. While the pollution 
outflow was more or less of 16 Tg in 2002 and 2004, the inflow in 2002 is 12% smaller than that 
in 2004. The four-year average pollution outflow is 17.7 Tg/year, of which about 25% (or 4.4 
Tg/year) can reach the west coast of North America. 

The interannual variation of pollution fluxes depends strongly on season and latitude, as 
shown in Figure 7 for the East Asia outflow. In the 30°-40°N segment, the springtime pollution 
flux in 2005 is about 180% higher than that in 2003, with that in 2002 and 2004 in between. The 
2005 summertime pollution flux is a factor of 2-3 higher than that of 2002 and 2004, and 30% 
higher than that of 2003. In latitudes north of 40°N, the springtime pollution flux in 2003 is 
about a factor of 2 in 40°-50°N and a factor of more than 3 higher in 50°-60°N than that in other 
years, due to record intense boreal forest fires. 

4.3. Estimated uncertainties of pollution fluxes 

The above-estimated pollution fluxes will inevitably be subject to uncertainties due to 
assumptions in implementing the algorithm. Major sources of uncertainties contributing to the 
pollution flux estimation integrated over latitudes and on annual average are listed in Table 3 
and discussed as follows. 
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One of first uncertainties would result from the MODIS -based estimate of anthropogenic 
aerosol optical depth. Kaufman et al. (2005b) estimated that the uncertainty in the derived 
anthropogenic AOD is 30%. Here we assume an uncertainty of 50% to reflect potentially larger 
uncertainty on regional scales. In this study seasonal and geographical variations of maritime 
fine-mode fraction have been accounted for. However, a robust validation or evaluation is not 
possible because of lack of observations of anthropogenic aerosol optical depth. 

The mass scattering efficiencies determined by aircraft and ship measurements are 
generally consistent ( Clarke et al., 2002) and their uncertainty comes from measurements of 
aerosol mass and scattering, which is estimated to be about 25% ( Quinn et al., 2004). The 
uncertainty for the single-scattering albedo is estimated to be 10%. These result in an estimate of 
overall uncertainty for the mass extinction efficiency to be 30%. For the fractional change in 
aerosol scattering as a function of RH, measurements from two nephelometers agreed within 
10% ( Carrico et al., 2003). For pollution aerosol, the hysteresis factor, a ratio of scattering 
coefficient on the monotonic to deliquescent curves of the hysteresis loop, is estimated to be 1 .3 
at RH = 60% ( Carrico et al., 2003). The use of average of monotonic and deliquescent growth 
curve would introduce an error of 15%. Therefore we estimate an overall uncertainty of 20% for 
f(RH). 


Polar-orbiting satellite sensors can only observe aerosols once daily and a combined use 
of the same MODIS sensor onboard the twin EOS satellites Terra (10:30AM local time) and 
Aqua (1:30PM local time) allows twice daily observations of aerosols in a region. Therefore a 
fundamental question to be addressed is: to what degree can flux estimates from MODIS/Terra 
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and/or MODIS/Aqua represent the diurnal average flux in all-sky conditions? While MODIS 
Terra and Aqua measurements on average well represent the daily aerosol optical depth on a 
global scale ( Kaufman et al., 2000, 2005c) and in East Asia ( Smirnov et al., 2002), it is not clear 
if MODIS-based estimates of pollution fluxes are a good representation of the daily average. 
Here we use GOCART model simulations to examine this issue by extracting the model 
calculated pollution fluxes around Terra/Aqua overpass times. It is found that the GOCART 
calculated fluxes at the Terra/ Aqua overpass times generally differ from the GOCART day and 
night averages by about 10%. These results give us some confidence that aerosol measurements 
from MODIS-Terra and/or MODIS-Aqua, can be used to quantify daily flux rates with a 
relatively small and definable uncertainty. 

The determination of transport heights of pollution aerosol plume has been used for dual 
purposes, i.e., selecting representative RH to account for humidification growth of aerosol and 
selecting representative wind speed for calculating mass fluxes. In this study, we have 
partitioned pollution AOT into 3 vertical segments of 0-1, 1-2, and above 2 km with equal 
fractions; a higher transport height (i.e., a larger fraction at higher altitudes) would result in 
larger aerosol mass flux because of stronger wind speed and drier airmass (hence greater aerosol 
mass for a given optical depth) at high altitudes in general. The sensitivity of pollution mass flux 
to the transport height should vary with region, depending on patterns of wind shear and 
atmospheric humidity. To examine such sensitivity, we assume pollution aerosol is evenly 
distributed in four vertical layers up to an altitude of 5 km and wind and humidity at 925, 850, 
700, and 600 hpa are used to calculate aerosol mass loading and fluxes in individual layers 
respectively. In comparison to the original assumption of vertical distribution, the pollution flux 
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increases by 20-80%, depending on season and latitude. By integrating over all latitudes, the East 
Asia outflow and North America inflow increases by 40-60% and 30-40%, respectively. 

The overall uncertainty in the calculated annual pollution aerosol fluxes for East Asia 
outflow and North America inflow is estimated as -80%, assuming that the individual 
uncertainties are independent. Note that the estimated pollution fluxes may have larger 
uncertainties in some sub-regions and seasons. While biomass burning smoke can come from 
both natural and manmade fires, it remains challenging to make an unambiguous distinction 
between them. Our assumption that all biomass burning aerosols are pollution would have biased 
our estimates of pollution fluxes in high latitudes in spring and summer where Eurasia boreal 
forest fires usually occur. On the other hand, the determination of transport heights and 
hygroscopic properties of pollution aerosol in this study is largely based on field experiments 
that have focused on industrial/urban pollution. This may have underestimated the flux of smoke 
aerosol because smoke aerosol can be raised to higher altitudes where wind speed is generally 
stronger, humidity is lower, and smoke aerosol usually has lower hygroscopicity (i.e., a larger 
mass loading for smoke than for industrial pollution for a given AOD). 

In this study we use the monthly mean dry aerosol mass derived from satellite retrievals 
in cloud-free conditions to calculate transport flux of pollution aerosol. Inherent in these 
estimates of pollutant transport is the assumption that the dry mass transport over the course of a 
month is the same, whether or not there are frontal cloud systems in the region. Frontal cloud 
systems extend over wide areas, preventing aerosol retrievals in those areas. These are not the 
same as broken cloud fields where satellite retrievals of aerosols are possible in the holes 
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between clouds. For example, because aerosol transport is associated with the WCB, and the 
WCB is associated with rising air in warm fronts that also consist of wide spread cloudiness, 
MODIS may miss important transport events every month. If dry aerosol mass in these events is 
significantly different from the monthly mean observed in cloud-free conditions, the estimated 
pollutant flux may have large uncertainties or biases. It has been suggested that a variety of 
compensating processes in and nearby clouds control the formation and removal of aerosols and 
hence determine their loading. While clouds can effectively remove aerosols from atmosphere 
through scavenging and rainout, they can also generate aerosol (e.g., in-cloud aqueous 
production of sulfate). We do not know how to quantify the differences in dry aerosol mass 
between cloud-free conditions and frontal cloud systems, and whether an overall bias exists in 
the estimated pollution fluxes. However, we acknowledge that the 80% uncertainty calculated 
from the known uncertainties may be a conservative number. 

A reduction of uncertainty in the estimated pollution aerosol fluxes can be achieved 
through integrating measurements from other A-Train sensors flying in formation with MODIS. 
For example, the new capabilities of passive remote sensing such as MODIS should be further 
enhanced by emerging lidar measurements of aerosol vertical distributions made by CALIPSO 
( Winker et al., 2003). GLAS data with limited temporal coverage has been used in this study to 
help determine aerosol transport heights. However, the availability of only 532 nm GLAS data 
makes it difficult to separate pollution aerosol from dust. Since June 2006, CALIPSO has been 
sensing the vertical distributions of aerosol backscatter and extinction in continuous mode, 
enhancing the temporal coverage of data. Furthermore the use of the two wavelengths and the 
depolarization capability at 532 nm of CALIOP can be used to better categorize aerosol types. In 
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the future, fully calibrated and validated CALIPSO aerosol extinction/backscatter profile data 
can be used to achieve a better quantification of aerosol intercontinental transport, in terms of 
seasonal, interannual, and geographical variations. A retrieval of fine and coarse-mode separated 
aerosol extinction profile with a fusion of CALIOP and MODIS/POLDER measurements 
( Kaufman et al., 2003a, b) would make a reliable estimate of vertical distributions of aerosol 
mass flux possible. 

The Polarization and Directionality of the Earth’s Reflectance (POLDER) on the 
Polarization and Anisotropy of Reflectance for Atmospheric Science coupled with Observations 
from a Lidar (PARASOL) is measuring directional and polarized radiance in nine spectral 
channels from 443 nm to 1020 nm. In addition to retrievals of AOD and size parameters, the 
analysis of the polarized angular scattered radiance provides additional information on the shape 
of particles (e.g., spherical vs non-spherical) ( Herman et al., 2005). This complements the 
MODIS spectral measurements which are not sensitive to particle shape and can be used to 
separate non-spherical dust from spherical maritime aerosols. POLDER and MODIS show an 
excellent agreement when particles are spherical ( Gerard et al., 2005). Therefore POLDER 
spherical coarse-mode aerosol optical depth over ocean can be used to separate dust from 
maritime aerosol in coarse-mode. The Ozone Monitoring Instrument (OMI) onboard Aura since 
June 2004 and its predecessor, the Total Ozone Mapping Spectrometer (TOMS) can separate 
UV-absorbing mineral dust or biomass burning smoke from scattering sulfate aerosol ( Herman et 
al., 1997). Such information could be used to determine if the plume is dominated by industrial 
pollution or smoke and hence to make a better characterization of aerosol growth with relative 
humidity. 
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5. Summary and Conclusion 

By taking advantage of MODIS high accuracy measurements of aerosol optical depth and 
fine-mode fraction over oceans, we have developed an observation-based approach to estimate 
the intercontinental transport of pollution aerosol over the North Pacific. For this method, 
measurements from other satellite sensors, such as AIRS and GLAS, and from field campaigns 
have also been used. We estimated that for 2002-2005 average, about 18 Tg/year of pollution 
aerosol is exported from Asia, of which about 25% arrives in the west coast of North America. 
The export and import pollution fluxes are largest in spring and smallest in summer, with the 
springtime and summertime flux accounting for about 40% and 10% of the annual integrated 
flux, respectively. For the period we have examined the strongest pollution export occurred in 
2003, due largely to record intense boreal forest fires in Eurasia. The overall uncertainty of 
pollution fluxes is estimated at 80%, due to uncertainties in separating pollution aerosol from 
dust and maritime aerosol, in deriving the aerosol mass loading and determining the transport 
heights of the pollution plume. A reduction of uncertainty can be achieved with a better 
characterization of pollution aerosol through integrating the PARASOL measurements of 
spherical and non-spherical aerosols, CALIPSO measurements of vertical distributions of 
pollution aerosol, and OMI aerosol indices. Simulations by the Goddard Ozone Chemistry 
Aerosol Radiation and Transport (GOCART) and Global Modeling Initiative (GMI) models 
agree quite well with the satellite-based estimates of annual and latitude-integrated fluxes, with 
larger model-satellite differences in latitudinal variations of fluxes. The satellite-based approach 
developed in this study can be adapted and extended to study the intercontinental transport of 
pollution in other regions. 
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While this study shows that a substantial amount of pollution aerosol from Asia can be 
transported to North America, its implications for climate change and air quality have not been 
assessed. Both modeling studies with reliable representations of major transport mechanisms at a 
wide range of scales and analyses of CALIPSO aerosol vertical distributions are necessary for 
gaining insights of downward transport of the height-elevated pollution to the atmospheric 
boundary layer and of potential influences of Asian pollution aerosol on surface air quality in 
North America. Further research is needed to assess how the intercontinental transport of 
pollution aerosol can influence the weather and climate over the North Pacific and how the 
resulting changes in atmospheric circulations would feedback on the transpacific transport of 
pollution and dust. Long-term trends in the transpacific pollution transport and its environmental 
impacts would be detected with an accumulation of satellite data in the future. 
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Table 1: Comparisons of MODIS-based estimate and GOCART and GMI simulations of 
East Asia outflow and North America inflow of pollution aerosol for 2004. 


Approach 

East Asia Outflow 
(Tg/year) 

North America 
Inflow 
(Tg/year) 

Inflow/Outflow 

(%) 

MODIS 

16.0 

4.2 

26 

GOCART 

15.0 

4.8 

32 

GMI 

14.4 

4.8 

34 


Table 2: Interannual variations of estimated outflow and inflow pollution fluxes 


Year 

East Asia Outflow 
(Tg/yr) 

North America Inflow 
(Tg/yr) 

Inflow/Outflow 

(%) 

2002 

15.8 

3.7 

23 

2003 

20.3 

5.7 

28 

2004 

16.0 

4.2 

26 

2005 

18.6 

4.1 

22 

2002-2005 

17.7 

4.4 

25 


Table 3: Estimated uncertainties in the flux calculation resulting from major 
assumptions in implementing the algorithm. The overall uncertainty is derived by 
assuming individual sources of uncertainty are independent. 


Major sources of uncertainties 

Estimated 

Uncertainty 

Derivation of pollution AOD 

50% 

Aerosol mass extinction efficiency (dry) MEE 

30% 

Humidification factor for the extinction f(RH) 

20% 

Using polar-orbiting measurements to represent 
diurnal average 

10% 

Transport heights 

50% 

Overall 

80% 
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Figure Captions: 


Figure 1: (a) Illustration of four boxes into which GLAS data are aggregated for deriving 
average profiles of aerosol backscatter and detection frequency during (b) September- 
November, 2003, and (c) February-March, 2004. 

Figure 2: Schematic of calculating pollution mass fluxes across meridional planes with a 
width of 10° in longitude in West Pacific (WP) and East Pacific (EP). The latitudinal 
integration over WP and EP is used to represent East Asia outflow and North America 
inflow, respectively. The background shows the distribution of pollution aerosol optical 
depth for MAM 2004. 

Figure 3: Seasonal variations of MODIS-derived East Asia outflow and North America 
inflow of pollution aerosol for 2004. 

Figure 4: Comparisons of seasonal variations of meridional integrated pollution flux 
across West Pacific and East Pacific. 

Figure 5: Comparisons of meridional variations of annual pollution flux across the West 
Pacific and East Pacific, respectively. 

Figure 6: Interannual variations of pollution aerosol optical depth for (a) MAM, and (b) 
JJA. 

Figure 7: Interannual variations of East Asia outflow pollution aerosol fluxes in different 
seasons and latitude segments. 
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Figure 2: Schematic of calculating pollution mass fluxes across meridional planes with a 
width of 10° in longitude in West Pacific (WP) and East Pacific (EP). The latitudinal 
integration over WP and EP is used to represent East Asia outflow and North America 
inflow, respectively. The background shows the distribution of pollution AOD for MAM 


2004. 



Figure 3: Seasonal variations of MODIS-derived East Asia outflow and North America 
inflow of pollution aerosol for 2004. 


44 





Transpacific Transport of Pollution Aerosol from Satellite 


Yu et al. 



Figure 4: Comparisons of seasonal variations of meridional integrated pollution flux 
across West Pacific and East Pacific. 



Figure 5: Comparisons of meridional variations of annual pollution flux across the West 
Pacific and East Pacific, respectively. 
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(a) 


MAM 2002 MAM 2003 



0.0 .05 0.1 .15 0.2 .25 0.3 0.4 0.5 1.0 


(b) 


JJA 2002 JJA 2003 



Figure 6: Interannual variations of pollution aerosol optical 
depth for (a) MAM and (b) JJA. 
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Figure 7: Interannual variations of East Asia outflow pollution aerosol fluxes in different 
seasons and latitude segments. 
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